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ABSTRACT 


This  paper  is  an  exposition  of  the,  GPSSR  system  to  support  management 
of  policy  and  execution  of  the  U.S.  Navy's  Enlisted  Personnel  Sea  Shore 
Rotation  Program.  Its  components  include  (1)  a  new  model  of  constrained 
network  goal  programming  type;  (2)  newly  developed  algorithms  for  use 
with  models  of  this  class;  (3)  computer  software  and  informatics 
developed  to  implement  these  algorithms,  plus  the  software  and  informatics 
for  other  modules  of  the  system  including  (4)  decision  support  tools 
for  report  generation  and  monitoring  capabilities. 
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1.0  INTRODUCTION 


This  report  presents  the  results  of  research  by  the  Center  for 
Cybernetic  Studies  to  provide  a  system  that  will  support  the  management  of 
policy  and  execution  of  the  Navy's  Enlisted  Personnel  Sea  Shore  Rotation 
Program.  This  system  consists  of  several  integrated  components  each  of 
which  represents  an  advance  in  the  present  state  of  modeling  and  computer¬ 
ized  algorithms.  These  components  include  (1)  a  new  model  of  constrained 
network  goal  programming  type;  (2)  newly  developed  algorithms  for  use  with 
models  of  this  class;  (3)  computer  software  and  informatics  developed  to 
implement  these  algorithms,  plus  the  software  and  informatics  for  other 
modules  of  the  system  including  (4)  decision  support  tools  for  report 
generation  and  monitoring  capabilities. 

The  system  was  developed  with  participation  by  the  staff  of  the  Navy 
Military  Personnel  Command  (NMPC)  1  and  the  Manpower  and  Career  Planning 
Research  Group  at  Carnegie-Mel Ion  University,  and  is  designed  to  meet 
present  Navy  requirements  for  both  planning  and  policy  evaluations. 


1  Special  thanks  are  due  to  CDR  E.  L.  Kainer  who  provided  guidance  and 
help  at  many  points  in  the  course  of  these  developments  and  to  ETCM 
William  F.  Hinkel  who  helped  in  collection  of  the  data,  formulation  of  the 
model,  and  interpretation  of  our  results.  LT.  Gareth  Habel  was  also 
helpful  at  critical  junctures  in  the  developments  covered  in  this  report. 
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Because  the  model  utilizes  a  goal  programming  approach,  we  refer  to  it  as 
the  Goal  Programming  Sea  Shore  Rotation  (GPSSR)  model. 


The  GPSSR  model  is  designed  for  use  in  planning  and  scheduling 
personnel  flows  and  for  evaluating  the  consequences  of  such  flows  relative 
to  Navy  goals  and  policies.  In  principle,  it  examines  all  possible 
personnel  flows  and  selects  the  ones  that  come  closest  to  meeting  all 
goals  while  honoring  the  specified  policy  and  operational  constraints.  It 
also  has  the  capability  of  evaluating  alternatives  in  policy  or  opera¬ 
tional  constraints  in  terms  of  their  effects  on  goal  achievements.  Thus, 
on  the  one  hand,  it  shows  what  is  required  to  do  the  best  possible  job 
under  the  given  constraints  and,  on  the  other  hand,  it  allows  the  explora¬ 
tion  of  alternatives  which  the  user  might  wish  to  consider.  By  providing 
a  consistent  basis  for  both  policy  and  operational  planning  through  its 
decision  support  tools,  GPSSR  also  provides  a  framework  for  policy 
execution  monitoring. 

+ 

Special  algorithms  developed  at  the  Center  for  Cybernetic  Studies 
which  have  now  been  incorporated  in  computer  software,  make  it  possible  to 
provide  the  above  capabilities  efficiently  and  effectively.  Solution 
times  of  a  minute  or  less  for  a  complete  detailing  community  have  already 
been  achieved  and  further  reductions  in  these  times  are  possible.  This 
provides  Navy  managers  with  a  capability  for  policy  analysis  and  planning 
at  various  desired  levels  of  detail  without  accompanying  delays:  a  Navy 
manager  may  try  a  variety  of  personnel  rotation  flow  alternatives  and  have 
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the  consequences  immediately  available  for  consideration.  The  report 
generation  and  graphic  display  capabilities  contained  in  the  system 
further  aid  the  manager's  decision  making. 


2.0  BACKGROUND 


The  management  of  Navy  enlisted  personnel  includes  the  continuous 

task  of  planning  for  and  executing  sea  shore  rotation  policies.  This 

management  task  is  described  in  the  Navy's  Enlisted  Transfer  Manual 

NAVPERS  15909C  Articles  3.0-3.01  as  follows: 

"The  system  for  the  planned  reassignment  of  personnel  among  the  vari¬ 
ous  types  of  duty  is  designed  to 

-  Promote  maximum  readiness  and  stability  both  afloat  and  ashore. 

-  Permit  equitable  opportunity  for  personnel  to  serve  in  duty  they 
consider  desirable. 

Rotation  among  sea,  shore,  and  overseas  activities  is  directly 
influenced  by  the  number  of  personnel  available  for  assignment, 
billets  authorized,  PCS  funds,  and  qualifications  of  the  individ¬ 
ual." 


Deciding  upon  rotation  policies  which  satisfy  a  variety  of 

oftentimes  conflicting  objectives  is  a  large  and  complex  problem  with  many 

different  dimensions.  Even  when  restricted  to  enlisted  personnel,  each  of 

more  than  250  detailing  communities  must  be  individually  considered,  and, 

for  effective  planning,  qualifications  like  the  following  are  involved: 

Paygrade 
Rati ng 

Subspecialty  (NEC)  Community 
Obligated  Service 
Contract  Group  or  LOS 
Individual  Starting  Point 
Prescribed  Sea  Tour 
Normal  Shore  Tour 


In  addition,  there  are  "exception  variables"  like  the  following: 

Special  Unit  or  Activity  Tour  Credit 

Unit  Deployment/Employment  Status 

Early  Release  Programs  Shipboard  Operational  Holds 

Permanent  Change  of  Station  (PCS)  Funding  Constraints 

Sex 

Voluntary  Shore  Duty  Curtailment 

Still  other  considerations  could  be  cited,  but  the  above  are  sufficient  to 
indicate  some  of  the  complexity  in  planning  sea/shore  rotations  and/or 
evaluating  policy  or  constraint  alternatives  for  their  sea/shore  rotation 
consequences . 

In  order  to  supplement  and/or  support  a  manually  operated  (or 
"stubby  pencil")  system,  several  unsuccessful  attempts  to  model  and 
computerize  the  sea/shore  rotation  process  were  undertaken.  The  first 
such  effort,  called  the  Dynamic  Flow  Model,  represented  an  attempt  at 
simulation  modeling  by  the  Navy  Personnel  Research  and  Development 
Center  (NPRDC)  in  the  early  1970's.  The  model  could  not  handle  a  sufficient 
number  of  the  essential  variables,  and  it  was  apoarent  that  efforts  to  extend 
and  enhance  these  capabilities  could  only  result  in  an  unwieldy  model. 

A  second  model  called  the  Discounted  Cash  Flow  (DCF)  Model  was 
developed  at  the  Office  of  the  Chief  of  Naval  Operations  (0P-01)  in  an 
effort  to  deal  with  "overall"  issues  of  policy.  At  this  level,  the  model 
failed  to  include  sufficient  detail  to  provide  any  real  insight  into  the 
rotation  problem. 
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A  third  effort  undertaken  in  considerable  detail  by  the  Center  for 
Naval  Analyses  (CNA)  was  completed  in  May  of  1979.  Called  by  a  variety  of 
names  —  CNA  Model,  Expanded  Sea/Shore  Rotation  Model,  ROTATIONMOD  —  this 
model  was  accepted  by  the  Navy  after  a  series  of  test  runs.  Partly  as  a 
result  of  changing  personnel  and  partly  as  a  result  of  subsequently 
discovered  deficiencies,  further  work  had  to  be  undertaken  in  order  to 
make  this  CNA  model  operational.  B-K  Dynamics,  Inc.,  was  retained  for 
this  work  and,  in  September  of  1982,  completed  a  user's  guide.  This  CNA 
model,  implemented  on  an  IBM/370  provc-d  to  be  slow,  expensive  and  confus¬ 
ing  to  use.  Quoting  from  [5]  1  ,  as  authored  by  B-K  Dynamics,  Inc., 

"The. .. system  is  expensive....  Please  keep  use  to  a 
minimum,  calling  up  the  model  only  when  a  course  of 
action  is  mapped  beforehand  and  a  computational 
strategy  designed.  This  will  cut  down  on  its  cost, 
which  could  be  surprisingly  high  when  the  system  is 
used  extensively." 

It  was  against  this  background  of  preceding  research  efforts  that 
the  work  on  GPSSR  was  undertaken.  More  than  a  system  for  effecting 

4 

Sea/Shore  rotation  was  intended.  By  agreement  with  the  Navy  and  the  CCS 
the  model  was  to  be  able  to  deal  with  rotation  scheduling  in  requisite 
detail  and  also  lend  itself  to  policy  evaluation  at  more  global  levels. 
It  was  also  to  provide  a  basis  for  improved  planning  of  "officer  based"  as 
well  as  "enlisted  based"  systems.  Finally  it  was  to  provide  a  possible 


1  Numbers  in  square  brackets  are  keyed  to  the  references  at  the  end  of 
thi s  report . 


approach  for  integrating  both  officer  and  enlisted  personnel  planning  to 
the  extent  that  this  might  be  feasible. 


Prior  experience  with  large  and  complex  personnel  planning  models 
made  it  clear  that  two  important  types  of  difficulties  were  to  be  anti¬ 
cipated  in  any  model  that  might  be  synthesized.  First,  a  variety  of 
conflicting  objectives  were  likely  to  be  encountered  so  that  some  way  was 
needed  for  dealing  with  the  difficulties  that  such  conflicts  can  cause  for 
most  types  of  mathematical  models.  "Goal  programming"  was  initially 
developed  by  A.  Charnes  and  W.  W.  Cooper  (in  collaboration  with  R. 
Ferguson  [2])  in  order  to  deal  with  such  conflicts  for  use  on  Navy  person¬ 
nel  problems.  Subsequently  extended  by  A.  Charnes  and  W.  W.  Cooper  (in 
collaboration  with  R.  Niehaus  [3]  and  [6])  it  also  has  the  capability  of 
showing  where  (and  in  what  amounts)  the  conflicts  are  causing  deviations 
from  prescribed  goals  and  policies. 

The  class  of  goal  programming  models  thus  provided  an  attractive 
basis  for  the  combinations  of  rotation  scheduling  and  policy  evaluation 
that  were  wanted.  This  was  one  reason  for  selecting  a  goal  programming 
approach  to  Sea/Shore  rotation.  Another  is  that  it  lends  itself  to  the 
kinds  of  extensions  that  might  subsequently  be  effected  to  "officer  based" 
as  well  as  "enlisted  based"  systems. 

A  second  class  of  difficulties  was  also  to  be  anticipated  in  the 
form  of  computational  algorithms  and  computer  codes  that  might  be  used  for 
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are  intended  only  to  show  the  sorts  of  reports  which  may  be  generated,  not 
the  results  from  a  real  detailing  community. 


Chart  1  was  obtained  by  downloading  data  from  the  network  optimizer 
to  an  IBM  PC,  then  graphing  this  data  with  the  LOTUS  1-2-3  program.  All  6 
duty  types  are  presented,  i.  e.  CONUS  shore  duty,  arduous  sea  duty,  over¬ 
sea  duty,  non-rotated  sea  duty,  neutral  duty,  and  oversea  prefered  duty. 
A  copy  of  the  instructions,  or  template,  for  the  IBM  PC  is  available  with 
the  GPSSR  system,  although  the  user  must  provide  a  copy  of  the  LOTUS  1-2-3 
program  in  order  to  use  this  template. 

Chart  2  accumulates  all  sea  and  shore  duty  so  the  user  may  see  the 
overall  Sea/Shore  picture.  This  chart  is  automatically  generated  by  our 
template  for  the  LOTUS  1-2-3  program  from  the  same  data  which  produced 
Chart  1. 
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in  Lovegren  [4],  the  current  version  has  achieved  another  order  of  magni¬ 
tude  increase  in  speed.  This  is  significant  in  its  own  right,  and  it  is 
also  indicative  of  the  progress  that  continues  to  be  made  in  reducing 
these  running  times.  Lovegren1 s  work  reduced  the  running  time  for  solving 
the  sea/shore  rotation  problem  from  24  hours  to  one  hour;  the  current 
version  runs  in  about  40  seconds  for  a  DC.  Furthemore,  the  previous 
version  did  not  take  into  account  "real  dollar"  costs,  as  does  the  current 
version.  To  distinguish  between  the  real  and  goal  costs,  the  new  code 
uses  an  approach  1  that  first  minimizes  deviations  from  stated  goals,  then 
achieves  this  result  at  the  lowest  possible  dollar  cost.  In  addition  the 
model  is  now  capable  of  keeping  track  of  different  kinds  of  dollars,  which 
can  be  important  when  funds  are  earmarked  and  non-transferabl e . 

4.5  REPORT  GENERATOR 

This  section  presents  a  summary  of  the  reports  which  may  be  obtained 
from  the  GPSSR  system.  These  output  modules  were  developed  concurrently 
with  the  modules  for  extracting  data,  so,  while  data  have  been  extracted 
from  the  Navy's  EMR  for  actual  detailing  communities,  these  report  modules 
were  tested  on  hypothetical  data,  and  the  charts  and  tables  presented  here 


1  This  corresponds  to  what  is  technically  called  a  "non-Archimedean" 
approach  as  described  in  detail  in  [I],  and  as  is  briefly  described  in 
Appendix  A. 
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ty,  costs  to  all  these  arcs,  so  that  it  makes  sense  for  the  program  to 
optimize  these  costs.  It  is  then  possible  to  obtain  the  set  of  flows 
which  minimizes  the  weighted  deviations  from  the  stated  goals  at  the  least 
possible  dollar  cost  while  remaining  within  the  constraints. 

At  the  present  time  a  file  has  already  been  written  with  a  set  of 
goal  and  dollar  costs.  This  has  been  done  so  potential  users  can  exper¬ 
iment  with  the  code  and  provide  possible  guidance  for  further  directions 
of  development.  Such  users  will  find  that  the  file  is  already  able  to 
provide  at  least  minimal  automatic  support  for  situations  in  which  the 
user  does  not  wish  to  supply  information  in  the  requisite  detail.  Users 
who  wish  to  do  so,  however,  can  insert  additional  information  about  costs, 
goals,  and  priorities  before  running  the  code.  Goals,  in  the  form  of 
manning  requirements  must  be  provided  to  the  model  at  this  time,  and  will 
be  used  to  write  "goal  arcs."  A  sample  file  is  available,  so  the  user  can 
see  the  proper  format  for  specifying  the  desired  billets.  A  description 
of  tjie  "goal  arcs"  used  to  represent  these  manning  requirements  is  in 
Lovegren  [4],  and  is  further  described  in  Appendix  A. 

4.4  NETWORK  OPTIMIZER  (VICNET) 

We  have  described  the  data  extraction,  computation  of  the  histor¬ 
ical  transition  rates,  and  the  generation  of  a  capacitated  network  with 
"valve  arcs,"  "bleeder  arcs,"  and  "goal  arcs,"  and  now  provide  a  brief 
description  of  the  network  optimizer.  As  compared  with  the  code  described 
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deviations  from  historical  rates,  but  only  with  a  penalty  cost.  It  is 
also  possible  to  maintain  rigid  (historical  rate)  constraints,  where 
these  are  known,  by  setting  5  to  0,  while  putting  prohibitive  penalties  on 
the  "bleeders."  A  more  detailed  description  is  deferred  to  Appendix  A. 

While  GPSSR  allows  a  great  deal  of  user  intervention,  it  is  designed 
so  that  it  does  not  place  heavy  burdens  on  the  user;  on  the  contrary,  very 
little  user  interaction  is  required.  The  user  need  only  call  the  appro¬ 
priate  optimization  routine  and  (optionally)  supply  flexibility  parame¬ 
ters.  The  rapid  solution  capabilities  of  the  optimization  algorithm  make 
it  feasible  to  explore  a  variety  of  alternatives  with  different  parame¬ 
ters.  Furthermore,  planned  enhancements  of  the  model's  user  interface 
will  largely  automate  this  process. 

4.3.2  ATTACHING  GOALS  AND  COSTS 

By  this  point,  a  network  has  been  generated  using  the  historical 
transition  rates  as  modified  by  the  user's  knowledge  and  experience.  A 
network  is  thus  obtained  with  arcs  which  describe  every  possible  transi¬ 
tion  from  one  paygrade,  length  of  service,  and  type  of  duty  to  some  other 
possible  combination  of  paygrade,  length  of  service,  and  type  of  duty.  On 
each  of  these  arcs,  we  have  also  imposed  upper  and  lower  bounds  which 
allow  flexibility  from  the  historical  proportions.  For  purposes  of  opti¬ 
mization,  it  is  then  necessary  to  attach  dollar  costs  and  goal,  or  priori- 
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projection  provides  an  estimate  of  the  flow  for  the  entire  period  covered 
by  the  model.  If  the  user  is  satisfied  with  such  a  quick  estimate,  and 
does  not  require  any  optimization,  it  is  possible  to  proceed  directly  to 
the  report  generator.  If,  however,  the  user  wishes  to  determine  the  flow 
of  personnel  which  will  "come  closest  to  meeting  goals  and  priorities"  at 
minimal  cost,  this  projection  will  then  be  embedded  in  a  constrained 
network. 

Part  of  the  flexibility  and  efficiency  of  GPSSR  comes  from  using 
this  projection  as  a  starting  point  for  developing  a  constrained  network. 
A  user-supplied  or  default  flexibility  parameter,  6  ,is  applied  to  the 
projected  flows  to  generate  upper  and  lower  bounds,  thus  allowing  flow  to 
occur  only  within  these  bounds  on  the  arcs  to  which  they  apply.  The 
resulting  network  is  a  pure  capacitated  network.  For  small  6  the 
constraints  are  satisfied  to  within  a  good  approximation,  even  where  the 
flows  within  the  indicated  bounds  do  not  satisfy  the  proportionality 
constraints  exactly.  Arcs^having  a  "window"  defined  by  such  upper  and 
lower  bounds  are  called  "valve  arcs."  The  flexibility  parameter  may  be 
varied  across  the  different  types  of  arcs,  so  that  windows  of  different 
sizes  can  be  generated  as  needed. 

The  approach,  as  described  to  this  point,  confines  the  model  to 
windows  determined  by  the  historical  rates  and  flexibility  parameter( s) . 
This  can  be  inappropriate  for  many  planning  and  evaluation  situations. 
Hence,  provision  is  made  for  the  addition  of  "bleeder"  arcs  to  permit 
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capacity  limit  which  the  flow  cannot  exceed.  The  introduction  of  these 
upper  and  lower  bounds  changes  the  model  from  a  pure  (or  uncapacited) 
network,  to  one  that  is  formally  characterized  as  having  a  capacitated 
network  structure. 


Ordinary  network  computer  codes  must  be  modified  to  deal  with 
networks  that  are  capacited.  GPSSR  must  also  handle  transition  conditions 
thct  involve  additional  "side"  constraints  so  that  still  further  exten¬ 
sions  of  these  network  codes  are  required.  We  have  avoided  the  use  of 
general  purpose  algorithms  for  networks  with  side  constraints — often 
called  a  "constrained  network" — because  these  algorithms  are  not  effi¬ 
cient  for  large  models  of  this  type.  For  models  as  large  as  ours  (for  a 
typical  DC,  a  network  with  several  thousand  nodes  and  tens  of  thousands  of 
arcs  is  generated,)  use  of  these  algorithms  requires  solution  times  which 
are  prohibitive. 


4.3.1  GENERATING  BOUNDS 


To  achieve  better  solution  times,  GPSSR  uses  a  new  algorithm  devel¬ 
oped  from  our  research  which  is  designed  to  approximate  this  "constrained 
network"  by  a  "pure  network"  which  is  also  capacited.  (A  more  mathemat¬ 
ical  description  of  these  types  of  networks  is  provided  in  Appendix  A.) 
This  is  done  in  two  steps:  First,  a  projection  routine  calculates  an 
exact  flow  on  each  arc  based  on  the  historical  transition  rates  generated 
by  the  previous  routine,  possibly  modified  by  an  informed  user.  This 
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4.2.3  USER  INTERACTION 

The  system  is  built  so  that  the  transition  rates  mentioned  above  can 
be  modified  by  the  user  for  those  cases  where  it  is  known  or  expected  that 
historical  transition  rates  will  not  reflect  the  actual  course  of  events. 
In  particular,  Enlisted  Community  Managers  generally  have  access  to  the 
planned  number  of  accessions  for  their  community,  a  number  which  may  be  at 
variance  with  the  historical  rates — e.g.  in  recent  years,  some  DCs  have 
experienced  significant  expansion  in  size.  For  these,  the  historical 
rates  of  accession  will  not  be  a  reasonable  indicator  of  the  actual 
accession  rates  observed.  When  this  occurs,  or  in  other  like  situations, 
the  user  can  input  the  planned  accessions,  overriding  the 
system-calculated  rates.  This  interaction  capability  is  currently  being 
upgraded  for  greater  "user  friendliness",  which  will  include  system 
supplied  prompts  and  menus  to  aid  users  in  their  choices. 

4.3  NETWORK  GENERATOR  * 

Having  described  the  modules  for  extracting  the  data  from  the  EMR, 
and  for  computing  the  smoothed  transition  rates  reflecting  the  historical 
proportions  of  promotion,  demotion,  etc.,  we  now  turn  to  the  third  module 
in  the  package,  the  network  generator.  This  module  introduces  upper  and 
lower  bounds  that  limit  the  personnel  flows  between  nodes  in  the 
network.  The  lower  bound  stipulates  a  minimal  amount  of  flow  that  must  be 
attained  on  the  arc  to  which  it  applies  while  the  upper  bound  provides  a 
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4.2.2  CALCULATING  THE  SMOOTHED  RATES 


After  determining  the  rank  and  type  duty  for  individual $,  and  how 
personnel  were  transferred,  this  module  takes  the  historical  transition 
totals  and  calculates  the  (Markov)  transition  rates  for  the  time  span 
covered  by  the  model.  This  is  accomplished  via  an  exponential  smoothing 
algorithm  which  uses  either  a  user-supplied  smoothing  factor,  or,  if  the 
user  prefers,  a  smoothing  factor,  a,  which  is  stored  in  the  computer.  The 
exponential  smoothing  algorithm  used  is  described  in  Appendix  B.  The 
existence  of  these  transition  rates,  as  reflected  in  the  proportionality 
constraints,  or  "side  constraints"  of  the  model,  would  normally  preclude 
solution  by  a  pure  network  program;  however,  by  relying  on  a  new  method  of 
approximating  these  constraints,  GPSSR  can  take  advantage  of  the  very  fast 
pure  network  codes  available  at  the  CCS.  The  new  method  is  explained  in 
detail  in  a  later  section. 
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satisfactory  separation  and  error  correction  have  been  achieved,  the 
remaining  modules  may  be  run  repeatedly  for  parameter  studies  without 
having  to  re-extract  this  data. 

4.2  CALCULATION  OF  TRANSITION  RATES 

Having  Ascribed  the  process  for  extracting  and  preparing  the  data, 
we  now  tur*i  to  the  second  of  the  five  GPSSR  modules,  the  transition  rate 
module,  which  computes  smoothed  Markov  rates  for  use  in  the  constraints 
for  the  network. 

4.2.1  OBTAINING  TRANSITION  TOTALS 

The  second  module  of  the  GPSSR  package  calculates  the  historical 
rates  of  accessions,  losses,  promotions,  demotions  and  rotations  in  the 
years  for  which  the  data  are  supplied.  This  is  done  by  examining  the 
extraction  from  the  EMR  for  two  successive  years,  finding  the  rank  anti¬ 
type  duty  for  individuals  in  the  DC  both  years,  and  thereby  determining 
how  many  individuals  were  promoted,  demoted,  rotated,  etc.  Individuals 
found  in  only  one  of  the  two  years  are  treated  as  accessions  or  losses  to 


4.0  GPSSR  SYSTEM  MODULES 


4.1  EXTRACTION  AND  SEPARATION  OF  DATA 

The  first  task  is  extraction  of  the  relevant  information  from  the 
Enlisted  Master  Record  (EMR).  The  raw  data  for  this  purpose  are  currently 
available  from  the  Center  for  Naval  Analyses  (CNA)  in  the  form  of  magnetic 
tapes.  Each  tape  contains  data  for  several  detailing  communities  (DCs) 
which  need  to  be  .separated  by  DC  for  use  in  this  model.  The  EMR  contains  a 
very  large  record  for  each  individual  from  which  only  a  few  data  fields 
are  needed.  After  these  fields  have  been  extracted,  a  DC-specific  file  is 
produced  containing  a  reduced  individual  record  for  each  member  in  the  DC. 

As  is  true  for  many  data  sources,  the  EMR  may  (and  generally  does) 
contain  some  errors  that  need  to  be  detected  and  corrected.  As  a  result, 
the  separation  programs  include  an  elaborate  structure  of  error-checking 
to  guarantee  "clean"  reduced  DC  files.  The  checking  is  accomplished,  in 
part,  through  use  of  the  many  fields  of  overlapping  information  found 
within  the  EMR.  Some  of  the  checking  cannot  be  done  automatically,  howev¬ 
er,  because  of  the  many  different  kinds  of  errors  potentially  to  be  found 
in  the  tapes,  and  so  the  programs  are  designed  to  enable  an  operator  to 
apply  his  or  her  own  knowledge  and  judgment  when  such  situations  arise. 
Even  so,  this  is  a  tedious  effort,  requiring  some  experience  with  the 
programs  as  well  as  knowledge  of  the  nature  of  the  data  in  the  EMR  tapes. 
The  separation  module  i s. Independent  of  the  other  modules;  hence,  once  a 
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then  computes  the  optimum  flows  on  this  network,  to  minimize  goal  costs  as 
well  as  dollar  costs.  Finally,  the  system  provides  a  report  generator  to 
display  various  aspects  of  this  optimal  solution  in  order  to  facilitate 
monitoring  and/or  redirection  of  these  efforts. 


1 
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Limitations  and  preferences  for  various  types  of  personnel  movement 
are  rendered  in  the  form  of  constraints  and  prescribed  goals  to  reflect 
given  rotation  policies.  Additional  constraints  include  the  transition 
rates  which  represent  the  historical  rates  of  promotion,  accession,  loss, 
etc.  of  personnel.  Two  different  kinds  of  goals  are  involved:  (1)  Those 
expressing  the  desire  to  fill  billets;  and  (2)  Those  expressing  the 
desire  to  rotate  personnel  in  accordance  with  Navy  priorities.  Goal  costs 
are  assigned  to  reflect  the  relative  importance  of  meeting  these  goals. 
Goals  are  derived  from  input  to  the  model  in  the  form  of  numbers  of  future 
personnel  authorizations,  or  proposed  changes  in  end  strength,  these 
changes  being  specified  as  numbers  or  as  percentages  of  current  staffing 
levels.  The  model  also  incorporates  the  real  dollar  costs  associated  vith 
the  Permanent  Change  of  Stations  (PCS)  involved.  Both  real  and  goal  costs 
are  reflected  in  the  minimizing  objective  as  noted  in  the  preceding  para¬ 
graph. 


The  GPSSR  system  consists  of  five  modules:  (1)  a  data  extraction 
component;  (2)  a  transition  rates  module;  (3)  a  network  generator;  (4)  a 
network  optimizer;  and  (5)  a  report  generator.  Mathematical  details  are 
supplied  in  Appendix  A  to  this  report.  The  operation  of  .GPSSR  may  be 
summarized  as  follows:  First,  the  model  extracts  raw  data  from  the 
Enlisted  Master  Record  (EMR).  Then,  after  providing  some  automated  data 
correction — as  well  as  the  facility  for  manual  data  adjustment — it  deter¬ 
mines  smoothed  historical  transition  rates.  When  these  transition  rates 
have  been  reviewed  and  respecified,  a  network  is  generated.  A  network  code 
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3.0  MODEL  AND  SYSTEM  OVERVIEW 


The  model  uses  a  goal  programming  network  form  for  representing  the 
flows  of  personnel  within  a  detailing  community  (DC)  over  time.  A  network, 
being  a  collection  of  nodes  and  arcs,  can  be  used  to  represent  states  and 
the  relations  of  flows  between  them.  In  GPSSR,  the  arcs  are  used  to 
represent  flows  of  personnel  between  the  nodes,  while  the  nodes  represent 
different  personnel  categories  and  status.  The  categories  are  defined  by 
those  qualifications  which  are  needed  to  capture  the  essence  of  the  rota¬ 
tion  problem,  and  it  is  the  number  of  these  qualifications  which  directly 
affects  the  problem  size.  Size  is  not  usually  a  problem  since  the  soft¬ 
ware  developed  by  the  CCS  is  presently  capable  of  handling  several  thou¬ 
sand  nodes  and  tens  of  thousands  of  arcs.  Some  understanding  of  the  model 
is  required,  however,  since  the  introduction  of  additional  parameters  can 
affect  the  size  of  the  problem  in  different  ways,  according  to  the  strate¬ 
gy  of  representation  used.  For  this  reason,  the  model  is  set  forth  in 
Appendix  A. 

Currently,  five  qualifications  are  used  to  define  the  nodes:  These 
are  an  individual's  (1)  paygrade;  (2)  time  on  tour;  (3)  length  of  service; 
(4)  type  of  duty;  and  (5)  the  specified  year  of  the  planning  horizon.  The 
objective  of  the  analytical  model  is  to  minimize  the  total  dollar  costs 
and  goal  costs,  subject  to  certain  constraints  and  network  relations. 
These  costs  and  constraints  are  discussed  below. 
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these  models.  Ordinary  goal  programming  computer  codes  would  not  be  up  to 
the  performances  required  in  these  applications.  Past  experience  with 
computer  codes  of  "network  varieties"  has  shown  that  these  types  of  codes 
can  now  accomodate  problems  of  huge  size  and  complexity,  provided  the 
problems  can  be  given  characterizations  that  lend  themselves  to  network 
representations.  Again,  A.  Charnes  and  W..  W.  Cooper  (in  collaboration 
with  R.  Niehaus  [3]  and  [6])  had  previous  experience  and  success  in  join¬ 
ing  network  and  goal  programming  models  in  a  single  goal  program¬ 
ming/network  representation  that  could  be  handled  by  available  network 
codes. 

In  the  present  case  (as  was  also  anticipated),  still  further  exten¬ 
sions  of  all  of  these  previous  developments  were  likelj^to  be  required. 
The  nature  of  these  extensions  are  described  in  the  sections  that  follow. 
For  clarity,  attention  is  confined  to  "enlisted  based"  Sea/Shore  Rotation 
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The  data  presented  by  Charts  1  and  2  may  also  be  obtained  in  tabular 
form.  However,  the  data  used  for  the  accompanying  table  is  not  the  same 
set  of  hypothetical  data  that  was  used  for  Charts  1  and  2,  since  the 
programs  to  generate  the  tables  and  charts  were  being  developed  in  paral¬ 
lel.  In  a  production  environment,  the  data  from  Table  1  would  be  down¬ 
loaded  from  a  mainframe  computer  to  an  IBM  PC,  or  some  other  personal 
computer,  and  input  to  our  template  for  the  LOTUS  1-2-3  program,  or  to 
some  other  program  with  similar  capabilities  to  produce  Charts  1  and  2. 

The  tabular  form  presents,  in  addition  to  the  information  in  the 
charts,  details  about  any  combination  of  scheduled  (i.e.  expected  under  an 
optimization  program)  promotions,  demotions,  accessions,  losses  and 
rotations  for  all  the  years  covered  by  the  model.  Table  1  presents  a 
sample  of  these  capabilities.  From  the  Table,  we  have  extracted  the  page 
presenting  CONUS  shore  duty,  arduous  sea  duty,  oversea  shore  duty,  and 
non-rotated  sea  duty  for  the  final  period  of  a  sample  run.  The  informa¬ 
tion  on  the  inventory  scheduled  by  the  optimizer,  the  user's  goals,  and 
the  deviations  from  those  goals  is  always  presented.  In  addition,  the 
user  requested  information  on  promotions,  losses,  and  accessions.  The 
total  movement  of  personnel  also  includes  demotions  and  rotations,  which 
were  not  requested  for  this  run  of  the  report  generator  but  which  can  also 
be  displayed.  As  a  result  of  this  flexibility,  the  vertical  columns  do 
not  sum  to  the  total  inventory  unless  all  categories  of  movement  are 
di splayed. 
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Table  2  contains  a  collage  of  the  larger  printout  from  which  the 
report  shown  in  the  previous  table  was  extracted.  This  printout  shows 
the  optimal  rotation  policy  broken  down  by  year,  type  of  rotation, 
paygrade,  time  on  tour,  and  length  of  service.  Several  sections  have  been 
pasted  together  to  present  a  better  view  than  was  possible  from  any  one 
section.  The  arc  numbers  and  names  indicate  the  different  sections  from 
which  they  were  extracted,  as  explained  below. 

At  the  top  of  the  table  is  shown,  as  the  problem  title,  the  name  of 
the  community  covered,  the  flexibility  option,  delta,  and  the  smoothing 
factor  alpha.  The  value  delta=0  shown  here  means  that  the  user  did  not 
use  the  flexibility  option,  and  the  smoothing  factor  alpha  =  0.2  was  used 
to  project  the  transitions  covered  by  the  exponential  smoothing  formula  of 
Appendix  B  in  this  case.  Finally,  the  total  goal  deviation  resulted  in 
penalties  of  15560  from  the  goal  deviation  penalties  used,  and  $280,020  is 
the  estimated  (best)  PCS  cost  associated  with  the  program  for  which  the 
details  in  Table  2  form  a  part. 
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Reading  from  left  to  right  the  column  headings  refer  to  the  follow 
ing:  ARCNUMBER 

This  is  the  number  of  the  arc  as  it  was  read  into  the  network  optimi¬ 
zer.  We  have  presented  a  selection  of  the  first  few  arcs,  and  three 
other  sections  taken  from  the  2000s  and  8000s.  The  first  few  arcs 
represent  initial  supply,  the  2000  arcs  represent  rotation  arcs 
(with  positive  dollar  cost),  the  first  set  of  8000  arcs  represent 
goal  arcs,  with  positive  penalty  costs,  and  the  second  set  of  8000 
arcs  represent  the  arcs  which  connect  the  goal  arcs  back  to  the 
beginning  of  the  network  to  form  a  complete  circuit. 

FROM  NODE 

These  are  the  source  nodes  from  which  each  arc  originates.  The  code 
tells  the  type  of  arc,  paygrade,  length  of  service,  etc. 

TO  NODE 

This  is  the  destination  of  the  arc.  1P3  02  means  (in  order)  year  1  of 
the  optimization,  promotion  arc  (P),  paygrade  E3,  length  of  service 
less  than  1  year  (0),  and  type  duty  2  (arduous  sea.) 

GOAL  COST 

Penalty  assigned  per  unit  flow  on  this  arc.  We  have  put  goal  costs 
in  this  formulation  only  on  failure  to  meet  desired  personnel  levels, 
with  the  -5  indicating  that  a  cost  of  5  units  is  assigned  to  falling 
below  the  requirement,  and  the  1  indicating  a  cost  of  1  unit  is 
assigned  to  exceeding  the  desired  level.  These  costs  are  not 
expected  necessarily  to  reflect  the  desires  of  actual  users.  Also, 
upper  and  lower  bounds  of  0  on  the  arc  with  -5  unit  cost  indicate 


29 


that  0  personnel  were  desired  for  this  category.  This  is  because, 
for  this  run,  no  goals  were  assigned,  so  the  program  used  0  for  all 
the  goals.  The  small  cost  of  exceeding  the  goals,  coupled  with  the 
lack  of  flexibility,  caused  personnel  to  be  scheduled  into  the  usual 
categories  anyway. 

DOLLAR  COST 

PCS  cost  per  person  assigned. 

UPPER  BOUND 

Maximum  flow  allowed  on  the  arc.  Since  no  flexibility  was  allowed 
(delta  =  0)  this  will  be  equal  to  the  lower  bound,  forcing  the  flow 
to  be  equal  to  the  set  upper  (or  lower)  bound,  except  on  certain 
"goal"  arcs,  where  violations  are  penalized  but  not  prevented. 

LOWER  BOUND 

The  minimum  flow  allowed  on  the  arc. 

ARC  FLOW 

Actual  flow  on  the  arc. 

ARC  COST  (G) 

Flow  multiplied  by  goal  cost. 

ARC  COST  ($) 

Flow  multiplied  by  dollar  cost. 

MARG  COST  (G) 

Marginal  cost,  i.e.  the  penalty  incurred  by  sending  one  more  person 
along  arc. 

MARG  COST  ($) 

Dollar  cost  incurred  by  one  more  person  along  arc. 
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5.0  MONITORING  FEATURES 


5.1  DISPLAYING  STATISTICS  FROM  THE  EMR 

Once  the  data  have  been  extracted  from  the  EMR,  the  University  of 
Texas  computer  system  provides  an  advanced  graphics  facility  which  makes 
it  possible  to  monitor  past  and  present  activities  and  consequences  of 
personnel  management  as  reflected  in  the  EMR  data.  This  is  an  important 
function  of  the  system,  because  the  size  and  complexity  of  personnel 
transfers,  as  well  as  the  existence  of  numerous  exceptions  often  masks  the 
real  situation  from  managers  trying  to  obtain  a  good  picture  with  only 
manual  methods  for  information  extraction  and  summarization  from  the 
data. For  example,  important  topics  like  how  much  of  an  existing  "rotation 
policy"  is  actually  being  implemented  in  view  of  the  exceptions  need  to  be 
addressed  regularly. 

As  a  start  toward  developing  desirable  monitoring  capabilities, 
GPSSR  currently  employs  the  statistical  package  (SAS)  to  calculate  and 
display  various  statistics  concerning  the  data.  The  package  can  produce  a 
graph  of  almost  any  combination  of  the  variables  found  in  the  EMR.  GPSSR 
provides  the  ability,  using  a  single  command,  to  generate  those  charts  and 
graphs  deemed  useful  for  policy  analysis.  The  following  examples  demon¬ 
strate  a  few  aspects  of  this  ability. 


32 


The  first  set  of  graphs,  represented  in  Figures  1  and  2,  is  obtained 
before  any  modelling  or  optimization  has  been  done.  In  principle,  the 
user  could  obtain  these  graphs  by  extracting  from  the  data  the  record  of 
every  individual  in  a  DC,  then  making  these  graphs  manually  by  plotting 
such  things  as  time  on  tour  vs.  type  duty  with  a  "stubby  pencil"  on  graph 
paper.  As  part  of  our  GPSSR  modeling  project,  however,  we  have  completely 
automated  the  process,  so  that,  with  a  single  command,  the  user  can  see 
these  results  for  purposes  of  monitoring  the  status  of  the  current  imple¬ 
mentation  of  Sea/Shore  rotation  policy  and  to  better  plan  future  rotation 
strategies.  For  example,  and  just  as  an  example,  we  have  chosen  to 
display  a  bar  chart,  showing  the  distribution  of  time  on  tour  for  two  of 
the  six  types  of  duty  in  DC  4000  at  Length  Of  Service  5-17  in  the  1982-1983 
time  frame.  The  two  types  of  duty  shown  are  (1)  CONUS  shore  duty,  in 
Figure  1;  and  (2)  Arduous  sea  duty,  in  Figure  2.  For  each  duty  type,  we 
display  a  histogram  showing  the  percentage  of  the  community  who  have  spent 
1,  2,  3,  4,  or  5+  years  assigned  to  that  type  duty  so  that  one  can  observe 
how  far  along  their  tour  most  of  the  community  lies,  and  where,  conse¬ 
quently,  an  extension  or  shortening  of  tour  length  will  have  the  most 


effect. 
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This  is  a  similar  histogram  for  oersonnel  on 
arduous  sea  duty  in  the  same  community.  By 
examining  the  distribution  of  personnel  with 
five  years  time  on  tour  the  user  can  see  the  effect 
of  a  small  extension. 
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Figure  2. 
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As  another  example  to  indicate  some  of  the  possibilities  of  such 
graphic  capabilities,  we  provide  the  following  plots  of  time  on  tour  vs. 
paygrade  in  Figures  3  and  4.  The  dotted  lines  show  the  distribution  of 
90%  of  the  community,  while  the  solid  line  shows  the  mean  over  all  lengths 
of  service.  Where  the  two  dotted  lines  diverge  very  markedly,  the  averages 
are  not  sufficiently  meaningful  for  drawing  firm  conclusions.  Converse¬ 
ly,  when  the  dotted  lines  lie  close  to  the  mean,  little  information  is 
lost  by  using  an  average  as  opposed  to  considering  all  the  observations 
separately.  Again,  we  show  these  plots  for  type  duties  1)  and  2),  CONUS 
shore  duty  and  arduous  sea  duty,  respectively. 
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This  is  a  graph  of  time  on  tour  vs.  paygrade 
for  personnel  on  CONUS  shore  duty.  The  dotted 
lines  enclose  90%  of  all  observations;  the 
sol  id  1 ine  is  the  mean. 
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Figure  3. 
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THE  GPSSR  NETWORK 


We  shall  shortly  present  a  simplified  numerical  illustration  but  first  we 
complete  our  network  interpretation  and  development  and  introduce  some  additional 
terminology  as  follows: 

In  GPSSR,  each  time  period's  network  is  divided  into  four  parts.  These  four  parts 
do  not  represent,  for  example,  the  four  seasons  of  a  year  or  the  four  weeks  of 
a  month,  but  are  just  logical  divisions  of  the  network.  The  first  section 
consists  of  the  promotion  spray  arcs.  These  arcs  handle  all  the  promotions, 
demotions,  losses,  and  accessions.  Type  of  duty  is  held  fixed  at  this  point. 

The  second  section  is  the  promotion  hose  arcs.  The  flow  through  a  node  is 
not  computed  by  the  network  optimizer,  only  flow  through  an  arc.  Every  node, 
representing  a  category  of  personnel,  is  then  connected  by  a  hose  arc  to  a 
node  representing  the  same  category,  so  the  flow  through  the  hose  arc  allows 
us  to  observe  the  flow  through  the  node,  which  is  equal  to  the  total  number  of 
personnel  in  that  category  during  the  time  period  in  question.  The  third 
section  consists  of  the  rotation  spray  arcs  which  handle  all  the  changes 
in  the  type  of  duty.  Finally  is  a  section  of  rotation  hose  arcs.  A  diagram  is 
given  in  Figure  (A-4). 

The  hose  arcs  have  not  been  altered  from  the  description  in  Lovegren  [4]. 

However,  the  promotion  spray  arc  (PR)  section  and  the  rotation  spray  arc  (RS) 
sections  of  that  network  have  been  modified  from  the  characterizations  used  by 
Lovegren.  It  is  impossible  to  show  an  entire  PR  section  or  RS  section  for  the 
actual  network.  However,  we  will  present  a  simplified  subsection  which  includes 
all  the  essential  features. 

We  first  present  a  simplified  PR  section  in  Figure  (A-5).  This  section 
is  replicated  for  every  combination  of  type  duty,  time  on  tour,  and  length 
of  service. 


Letting  _T ,jj  *  P  ±  6)  ,  ( A- 11)  is  seen  to  be  like  (A-5)  and  we  may 
employ  the  power  of  our  specialized  network  optimizer. 

NON-ARCHIMEDEAN  NETWORK  OPTIMIZER 

The  cost  vector  £  in  ( A- 11)  represents  a  set  of  penalties  for  failing  to 
achieve  goals.  They  are  artifacts  of  the  model  and  not  actual  dollar  costs. 

In  some  scenarios,  it  is  imperative  that  goals  be  met  as  closely  as  possible 
regardless  of  cost,  but  if  alternate  solutions  exist  which  meet  the  goals  e- 
qually  well,  then  the  solution  which  minimized  real  dollar  costs  should  be  chosen. 
This  is  achieved  by  using  a  non-Archimedean  (or  non-standard)  version  of  the  net¬ 
work  optimizer,  which  solves 


min  cjx  +  e£jx 
s.t.  Nx  =  a 

]_<£<£ 

where  e  is  a  non-standard  infinites  imal  In  this  formulation, 
emptively"  minimized  which  means  that  c^x  is  considered  only  in 
not  alter  optimal  values  of  cjx..  With  this  formulation  the 
mizer  will  achieve  stated  goals  as  closely  as  possible  and  then 
costs 


( A- 12 ) 


£^x  is  pre- 
a  way  that  will 
network  opti- 
minimize  dollar 


1 

Alternatively,  if  these  are  two  sets  of  goals,  one  pre-emptively  important, 
then  two  sets  of  penalty  costs  could  be  used. 


Algorithms  to  solve  problems  of  the  form  (A-9)  exactly  are  two  orders  of 
magnitude  slower  than  algorithms  to  solve  (A-5).  For  GPSSR,  we  assume  that  the 
rates  of  promotion,  loss,  etc.  will  be  similar  to  the  historical  rates,  and  these 
are  our  side  constraints.  Thus,  the  network  constraint  on  the  flow  through  a 
node  j  is  simply  the  previously  discussed  condition  that  flow  in  equals  flow 
out.  The  proportionality  constraints  associated  with  P  put  an  additional  re¬ 
quirement  on  the  flows  out  of  the  node--e.g.  they  must  be  proportional  to  the 
total  flows  through  the  node. 

It  is  reasonable  to  solve  (A-9)  approximately.  In  fact,  since  we  do  not 
expect  historical  rates  to  be  followed  exactly,  a  more  realistic  version  of 
(A-9)  is 

min  cTx  (A-10) 

s.t.  Nx  =  a 

-6  <  Px  -  b  <  o 

where  the  components  of  6  represent  the  maximum  and  the  components  of  -6 
represent  the  minimum  admissible  deviations  from  the  corresponding  components 
of  b. 


By  assigning  a  proportionality  constraint  to  every  arc,  P  becomes  inver¬ 
tible.  Then  we  may  write 

min  cTx  (A-ll) 

_x 

s.t.  Nx^  =  a 

P_1(b  -6)  ^  x  P-1(b  +  5) 
since  P~*is  non-negative. 
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FIGURE  A-3 


NETWORKS  WITH  SIDE  CONSTRAINTS 

Let  N  be  a  node  incidence  matrix,  P  an  arbitrary  matrix,  and  consider  the 
problem 

min  c  x  (A-9) 

x 

s.t.  Nx  =  _a 
P  x  =  b 
J_  <  x.  <  u. 

This  is  a  network  with  side  constraints,  the  side  constraints  being  Px  =  b, 
where  P  is  a  matrix  of  coefficients  and  b  a  vector  of  additional  conditions.  In 
addition  the  1  and  u  are  vectors  that  impose  lower  and  upper  bounds  on  the- pos¬ 
sible  choices  of  x. 
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Minimizing  this  function  gets  the  x. .  "as  close  as  possible"  to  the  g... 

ij  3ij 

The  objective  function  (A-7),  however,  does  not  appear  to  be  in  the  form 
required  by  (A-5).  A  transformation  due  to  Charnes,  Cooper,  et  al .  [l] 
brings  us  back  to  (A-5): 


mi 


inZZcij(5ij +  6ij> 

i  J 


(A-8) 


s.t.  x.  6.  .  +  6~. .  =  g.  ., 
ij  ij  ij 


Note  that  the  <5..,  represent  deviations  above  and  below  the  goal  g..  resulting 
from  the  flow  value  assigned  to  x^..  These  deviations  are  accorded  penalties 
or  "costs"  c.j  per  unit  in  the  functional  being  minimized. 

The  network  representation  of  (A-8)  in  terms  of  goal  arcs  is  shown  in 
Figure  (A-2),  and  an  alternative  but  equivalent  representation  in  Figure  (A-3). 


FIGURE  A-2 
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than  a  general  purpose  LP  package  such  as  MPSX  when  applied  to  an  optimization 
of  the  form  (A-5). 

REPRESENTATION  OF  NETWORKS 

Given  a  problem  in  the  form  (A-5)  it  is  easiest  to  visualize  the  problem 
by  drawing  it  as  a  network.  An  arc  of  the  network  is  represented  as  in 
Figure  (A-l),  where  c..  is  the  unit  cost  associated  with  the  arc,  1. .  is  the 

1  J  *  J 

lower  bound,  or  minimum  flow  required  on  the  arc,  and  u. .  is  the  upper  bound, 

*  J 

or  maximum  flow  which  may  be  allowed.  (It  is  impossible  to  display  the  en¬ 
tire  GPSSR  network,  as  it  contains  over  10,000  nodes  and  20,000  arcs.  Enough 
simplified  subsections  are  presented  below  to  give  a  good  picture  of  the  en¬ 
tire  network.) 

FIGURE  A-l 


GOAL  PROGRAMMING 

The  GPSSR  program  tries  to  meet  personnel  goals  as  closely  as  possible. 
This  implies  that  the  objective  function  contains  terms  of  the  form 

f 1  •  9iji 

ij  i  j 

where  g^  is  the  goal  and  the  vertical  strokes  represent  an  absolute  value  of 

the  difference  between  g..  and  x... 

ij 
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each  row  associated  with  a  node  for  which  non-zero  entries  appear.  Each  column 
represents  an  arc  with  the  ±1  values  indicating  the  nodes  on  which  it  is  incident. 
That  is,  since  the  column  has  a  +1  in  row  j,  say,  and  a  -1  in  row  i,  it  may  be 
graphically  represented  as  an  arc  from  node  i  to  node  j. 

Additional  constraints  of  the  following  type  may  be  added  to  form  a 
(pure)  capacitated  network  problem: 


1  . .  <  x . .  <  u . . 

U  iJ  iJ 


(A-4) 


The  complete  problem  is  then 


min  cTx 


(A-5) 


s.t. 


Nx  =  a 
1  <  x  <  u 


Note  that  (A-5)  looks  similar  to  the  general  linear  programming  (LP) 
problem 


min  cTx 


(A-6) 


s.t. 


Ax  =  jj 
x  >  0 


where  A  is  an  arbitrary  matrix.  However,  the  algorithms  to  solve  (A-6)  require 
some  two  orders  of  magnitude  more  computations  than  the  algorithms  to  solve 
(A-5).  In  particular,  the  Center  for  Cybernetic  Studies  has  developed  one 
of  the  most  efficient  network  optimizers  available,  a  specialized  package 
which  can  solve  (A-5)  but  not  (A-6),  and  is  two  orders  of  magnitude  faster 
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MATHEMATICAL  DESCRIPTION  OF  THE  MODEL 

Mach  of  the  following  is  abridged  from  Lovegren  (for  a  fuller  explanation 

see  [4]).  A  network  may  be  visualized  as  a  collection  of  nodes  S  =  {1,2 . n}, 

and  between  these  nodes  a  set  of  arcs.  Along  each  arc  is  a  flow  x.^. ,  the  flow 
from  node  i  to  node  j.  If  x..  <  0,  this  represents  a  flow  of  |x..|  from  node 

•  J  ’  J 

j  to  node  i.  Using  c..  to  represent  the  cost  per  unit  flow  from  node  i  to  node 

J 

j,  the  pure  network  optimization  problem  is  then 

(A-i) 

ij  i  J 

subject  to  the  network  constraints  "lo  que  entra  sale"  or  "what  comes  in  goes 
out,"  i.e.. 


Dxik  -  Exkj =  ak  <A-2) 

i-  J 

This  says  that,  at  each  node  k,  the  total  flow  going  into  the  node  minus  the 
total  flow  going  out  of  the  node  is  equal  to  the  net  inflow  or  outflow  at  that 
node.  In  matrix  notation,  (A-l),  (A-2)  can  be  written 

min  £Tx  (A-3) 

s.t.  Nx  3  a 

Components  of  the  c  vector  represent  the  cost  per  unit  flow  on  each  arc  and  the 
component  of  the  21  vector  represent  these  flows  (from  node  i  to  node  j).  Since 
every  arc  must  go  between  two  nodes--into  one  and  out  of  the  other--each  column 
of  the  N  matrix  has  precisely  two  non-zero  entries:  +1  and  -1.  All  other  entries 
in  each  column  are  zero  except  for  these  ±1  values  which  are  incident  on  nodes  i 
and  j,  respectively.  The  N  matrix  is  called  a  node-incidence  matrix.  In  fact, 
any  matrix  with  this  property  may  be  considered  a  node  incidence  matrix,  with 
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7.0  SUMMARY  AND  CONCLUSIONS 

7.1  SUMMARY 

The  GPSSR  system  is  a  sophisticated  Management  Information/  Deci¬ 
sion  Support  System,  produced  by  the  Center  for  Cybernetic  Studies  at  the 
University  of  Texas  for  the  U.  S.  Navy.  The  system  handles  possibly 

contradictory  information  by  optimizing,  via  goal -programmi ng ,  over  suit¬ 
able  grils,  using  a  capacitated  network  model  structure  with  computation 
orders  of  magnitude  faster  than  that  of  previous  Sea/Shore  rotation 

models.  The  system  contains  a  monitoring  capability  which  provides  a 
manager  with  previously  unavailable  information  about  the  Sea/Shore  rota¬ 
tion  policy  actually  implemented. 

7.2  CONCLUSIONS  AND  DIRECTIONS  OF  FURTHER  WORK 

This  system  will  be  a  useful  tool  for  Navy  managers  and  planners. 
It  is  also  general  enough  to  be  applied  in  solving  an  array  of  problems 

other  than  the  sea/shore  rotation  problem.  It  can  be  used  to  solve  any 

problem — including  optimization  problems — with  elements  involving  goals 
and  flows,  capacities  and  costs.  Its  goal  programming  features  permit 
identification  and  analyses  of  deviations  from  goals  caused  by  one  or  more 
of  these  elements  such  as  might  be  involved  in  officer  or  enlisted  based 
problems  and  the  planning  of  optimal  force  structures.  Finally,  the  dual 
evaluators  are  available  for  exploitation  in  policy  analyses  and  evalu¬ 
ations  such  as  are  likely  to  be  present  in  allocation  policy  problems 
associated  with  manpower  planning. 
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6.2  TECHNICAL  ENHANCEMENTS/REFINEMENTS 


•  Developing  an  "intelligent"  user  interface  for  the  GPSSR.  One  of 
the  main  goals  of  this  effort  is  to  provide  powerful  interactive 
capabilities,  so  that  a  decision  maker  need  be  neither  a  computer 
expert  nor  an  operations  research  expert  in  order  to  use  the 
system.  Using  normal  Navy  language,  the  user  should  be  able  to 
explain  his  problem  to  the  system,  which  will  automatically  call 
the  appropriate  programs,  prompt  the  user  for  specification  of 
parameters  and  directives,  and  produce  the  desired  output. 
Natural  language  processing  in  all  detail  is  more  ambitious  than  we 
expect  to  achieve,  but  we  do  intend  to  push  very  far  in  that  direc¬ 
tion  . 

•  Producing  more  summary  reports  as  derived  from  the  global  output 
f i le.One  such  report  should  aggregate  the  costs  resulting  from  the 
personnel  movement  in  the  network.  Currently,  the  aggregation  is 
by  paygrade  and  type  duty,  calculated  separately  for  the  different 
types  of  costs. 

•  Enhancing  the  guality  of  the  input  data.  We  hope  to  improve  our 
understanding  of  how  to  handle  some  of  the  inputs  which  have  not, 
as  yet,  been  thoroughly  checked.  The  quality  of  cost  information, 
for  example,  must  be  improved. 


6.0  WORK  IN  PROGRESS 


Parallel  efforts  are  also  currently  under  way  in  the  CCS  to  improve 
the  GPSSR  capability  and  performance.  A  brief  description  of  some  of  these 
follows  :  - 


6.1  DEVELOPMENTS  IN  THE  THEORY 


Introducing  an  alternative  goal  concept,  that  of  goal  "length  of 
tours . 11  To  that  end,  the  time  on  tour  has  been  introduced  as 
another  node  dimension,  enabling  the  model  to  calculate  penalties 
based  on  deviations  of  desired  lengths  of  tour.  These  penalties 
are  then  added  in  with  the  other  goal  costs,  representing  devi¬ 
ations  from  the  planned  billets,  which  were  already  in  the  model. 
Early  rotations,  which  might  be  of  concern  to  the  DC  personnel 
management,  are  penalized.  Likewise,  late  rotations,  which  might 
cause  individuals  to  quit,  are  also  penalized. 

Studying  the  effects  of  the  non-Archimedean  optimization  on  the 
rate  of  change  in  the  DC  strength. As  explained  in  3.4,  the  minimi¬ 
zation  is  taken  first  on  the  goal  costs  and  only  then  over  the  real 
costs.  If  the  end  strength  goals  are  somewhat  higher  than  the  start 
inventory,  the  model,  given  only  the  dollar  costs  for  maintaining 
personnel,  and  only  goals  for  strength  in  the  final  year,  will  try 
to  access  people  as  late  as  possible  to  avoid  the  costs  of  carrying 
them  along  the  network.  The  computed  solution  may  then  suggest 
abrupt  changes  in  manning  for  the  DC,  all  taking  place  in  the  last 
year  under  consideration.  However,  the  introduction  of  interme¬ 
diate  goals  via  "valve  arcs"  and  "bleeder"  arcs,  will  cause  the 
model  to  provide  for  gradual  changes  and  smooth-out  possible 
saw-like  jumps  in  the  personnel  curve. 

Considering  different  scenarios  and  objectives  regarding  the  male 
and  female  personnel  in  certain  DCs.  Special  attention  needs  to  be 
devoted  in  the  modeling  process  to  address  problems  resulting  from 
legal  constraints  and  the  lack  of  available  positions  at  sea  which 
are  adequate  for  women  (e.g.  older  ships  must  be  modified  to  accom¬ 
modate  female  personnel).  This  situation  creates  imbalances  in 
the  rotation  policies  applicable  to  different  sexes. 
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PERCENTAGE  PROMOTED  VS  PAYGRADE  FOR  TYPE  DUTY  2 

DETAILING  COMMUNITY  4000 


This  is  a  graph  representing  a  hypothetical 
community  showing  the  rate  at  which  personnel 
are  promoted  as  a  function  of  paygrade.  The 
data  do  not  represent  a  real  or  necessarily 
representative  community.  The  promotion  rates 
are  averaged  over  all  lengths  of  service  between 
5  and  18,  and  the  solid  line  shows  the  mean, 
while  the  dotted  lines  enclose  90%  of  all 
observations 


PAYGRADE 

LENGTH  OF  SERVICE  BETWEEN  5  -  18 
Figure  5- 


5.2  ADDITIONAL  MONITORING  CAPABILITIES 


Not  immediately  available  from  the  EMR  are  rates  of  accession  and 
loss,  and  length  of  tour  as  opposed  to  time  on  tour.  In  order  to  obtain 
these  rates,  we  had  to  compare  two  years  of  the  EMR.  Note  that  accession 
and  loss,  for  our  purposes,  refer  to  a  single  community.  People  who 
transfer  from  one  DC  to  another  are  considered  an  accession  to  their  new 
community  and  a  loss  to  their  old  community.  For  purposes  of  filling  a 
given  community's  billets,  this  should  not  be  an  unreasonable  definition. 

These  data,  as  well  as  promotion  and  demotion  rates  are  available  in 
a  readable  file,  and  plans  exist  for  a  report  generator  that  will  make 
them  even  more  accessible.  In  addition,  a  graph  package  is  planned  that 
will  present  the  data  in  a  form  similar  to  the  example  shown.  The  example 
was  prepared  using  the  SAS  package.  However,  as  part  of  the  continuing 
effort  to  develop  an  intelligent  user  interface,  a  more  user-friendly  plot 
interface  is  planned  which  will  be  much  easier  to  access  than  the  SAS  plot 
package . 

As  a  further  example  of  the  kinds  of  GPSSR  graphs  that  can  be  gener¬ 
ated  at  this  point,  the  chart  in  Figure  5,  based  on  hypothetical  data,  is 
presented  to  show  the  proportion  of  personnel  promoted  while  on  type  duty 
2  as  a  function  of  paygrade.  Here  the  dotted  lines  show  the  range  of  90% 
of  the  data,  averaged  over  length  of  service  and  time  on  tour. 
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Figure  4. 


FIGURE  A-4 


CONTINUATION  BLEEDER  ARC 


FIGURE  A-5 


This  is  still  not  a  complete  diagram  of  a  PR  subsection,  as  the  accession 
and  loss  arcs  have  been  omitted  for  purposes  of  simplifying  the  representation. 
The  section  illustrates  how  valve  arcs  and  bleeder  arcs  are  used  to  account 
for  historical  rates  of  promotion,  demotion,  accession  and  loss.  Flowing  into 
this  section  are  the  RS  hose  arcs,  and  the  flow  continues  with  the  PR  hose  arcs 
neither  of  which  are  shown.  What  is  shown  are  the  valve  and  bleeder  arcs  for 
promotion  and  demotion  (accession  and  loss  being  handled  similarly).  The 
valve  arcs  allow  personnel  to  move  through  the  network  at  historical  rates 
(±6)  with  no  penalty  costs.  If  historical  rates  are  to  be  violated  by  more  than 
a  specified  percentage,  however,  the  bleeder  arcs  assign  appropriate  pen¬ 
alties  which  will  be  incurred  in  the  violation.  These  penalty  rates  are  used 
to  discourage  violations  and,  indeed,  increasing  deviations  can  be  penalized 
at  increasing  values— although  this  is  not  illustrated  in  the  diagram. 

Iri  the  actual  diagram  for  a  hypothetical  community,  45  people  would  have 
been  promoted  from  E3  to  E4  based  on  historical  rates.  Thus,  the  segment 
drawn  allows  40-50  persons  to  be  promoted  with  no  penalty,  but  imposes  a  cost 
of  five  units  per  person  for  deviations  above  or  below  this  range. 

We  next  show  a  simplified  version  of  the  RS  spray  arcs  for  the  last 
period  of  the  problem,  along  with  the  goal  arcs  in  Figure  (A-6).  We  have 
shown  only  a  single  paygrade  and  length  of  service.  Also,  coming  into  the 
section  are  the  promotion  hose  arcs,  which  are  not  shown.  We  illustrate  a 
situation  with  a  maximum  tour  length  of  three  periods,  and  only  three  types 
of  duty.  Even  so,  this  subsection  has  18  nodes  for  the  RS  part,  and  27  RS 
spray  arcs,  only  9  of  which  have  been  drawn. 

We  assume  that  policy  is  to  keep  personnel  in  duty  type  three  for  three 
years,  then  transfer  them,  if  possible,  to  duty  type  two.  If  this  is  not 
possible,  the  second  choice  is  a  transfer  to  duty  type  one,  with  the  last 


FIGURE  A- 6 


choice  being  for  personnel  to  remain  an  additional  period  in  duty  type  three. 
Penalties  are  therefore  assigned  to  premature  rotations,  to  late  rotations,  and 
to  rotations  to  duty  type  one  rather  than  duty  type  two,  so  that  three  years  of 
duty  type  three  would  give  the  minimum  amount  of  penalty  charges.  The  rela¬ 
tive  weights  are  purely  hypothetical,  chosen  for  illustration  only. 

We  have  assumed  that  time  on  tour  is  not  relevant  to  meeting  the  staffing 
levels  indicated  by  the  goal  arcs  in  the  diagram,  so  converge  arcs  sum  over¬ 
all  times  on  tour  to  the  converge  nodes  which  are  then  connected  by  the  ter¬ 
minal  nodes  by  the  goal  arcs.  Goal  costs  and  upper  and  lower  bounds  exist 
on  all  the  goal  arcs,  but  are  only  written  on  the  first  set  of  goal  arcs  in 
the  diagram. 

A  Numerical  Example 

In  this  section,  we  shall  present,  a  simple  illustrative  example. 

For  this,  let  us  assume  that  the  community  is  stable,  i.e.  that  the 
number  of  personnel  lost  equals  (approximately)  the  number  gained 
through  promotions,  demotions,  attritions,  etc.  Let  us  then  concentrate 
on  the  rotation  policy  goals  vs.  staffing  goals  for  a  single  paygrade 
i.e,  we  shall  fill  in  the  data  for  the  arcs  for  Type  Dutys  1  and  2  in 
Figure  A-6,  and  solve  the  resulting  network..  Recall  from  page  A10,  that 
the  first  choice  for  Outy  Type  3  was  to  rotate  to  Duty  Type  2,  and  the 
second  choice  was  to  rotate  to  Duty  Type  1.  For  our  example,  suppose 
that  the  first  choice  for  Duty  Type  2  is  to  rotate  to  Duty  Type  1,  with 
second  choice  being  to  rotate  to  Duty  Type  3,  and  that  the  first  choice 
for  Outy  Type  1  is  to  rotate  to  Duty  Type  3,  with  Duty  Type  2  the  second 
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choice.  These  choice  preferences  are  summarized  in  Table  A-l. 


Type  duty  To 

1 

2 

3 

From 

1 

3rd 

2nd 

1st 

2 

1st 

3rd 

2nd 

3 

2nd 

1st 

3rd 

Table  A-l  Goals  for  type  of  rotation 

The  above  table  gives  the  preferred  rotation  sequence.  The  next 
consideration  is  tour  length.  We  assume  that  policy  is  to  have 
personnel  serve  1  full  period  of  Type  1  Duty,  and  two  full  periods  of 
Type  2  anG  3  Duty  before  transfer,  as  summarized  in  Table  A-2. 

Type  'uty  1  2  3 

Tour  Lenith  122 

Table  A-2  Goals  for  tour  lengths 


Finally,  we  assume  that  the  iverage  PCS  cost  is  $5000  from  Duty 
Type  1  to  Duty  Type  2,  $9000  from  Du  v  Type  2  to  Duty  Type  3,  and 
$20,000  from  Duty  Type  1  to  Duty  Type  3,  as  summarized  in  Table  A-3. 


Type  duty  To 

1 

2 

3 

From 

1 

SO 

$5000 

$20,000 

2 

S5000 

$0 

$9000 

3 

$20,000 

$9000 

so 

Table  A-3  Average  PCS  Costs 


i 
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TYPE  DUTY  2  TYPE  DUTY  3 


J 


These  data  are  hypothetical,  of  course,  and  greatly  simplified  for 
ease  of  understanding,  to  only  three  types  of  duty  and  three  periods  for 
tour  length.  The  part  of  these  data  for  Type  Duty  3  are  drawn  in  the 
network  fragment  of  Figure  A-6,  on  page  All.  An  extract  of  a  single  arc 
of  this  figure  is  presented  as  Figure  A-7.  On  this  arc  from  Type  Duty 
3,  time  on  tour  2  to  Type  Duty  2,  time  on  tour  1,  we  have  indicated  the 
2  units  of  goal  cost  and  also  the  $9000  PCS  cost  as  represented  in  the 
above  tables.  See  the  circled  numbers  on  this  arc.  The  parenthesized 
values,  (0,~)  mean  that  all  flows  are  in  the  direction  indicated  by  the 
arrow  and  there  is  no  upper  bound  imposed  on  these  flows. 


TIME  ON  TOUR  2 


I 

i  •  •  #  • 


TIME  ON  TOUR  1 

•  •  •  • 

•  •  •  • 


FIGURE  A-7 

Legend:  The  dots  indicate  the  presence  of  other  node-arc  incidences  and 
penalties  which  are  shown  in  detail  in  Figure  A-6. 
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Again  we  emphasize  that  we  are  trying  to  restrict  this  discussion 
to  simple  versions  of  a  complex  problem— -whi le  recalling  that  our  model 
and  algorithm  with  associated  software  can  handle  very  large  problems 
with  extremely  fast  solution  times.  Since  drawing  the  entire  network  of 
24  nodes  and  42  arcs  would  only  complicate  our  discussion  for  this 
example,  the  arcs  for  Types  Duty  1  and  2  are  merely  indicated,  rather 
than  drawn  in  full. 

Referring  to  Figure  A-6,  the  set  of  arcs  on  the  extreme  right  are 
GOAL  arcs;  the  arcs  adjacent  to  the  GOAL  arcs  are  the  CONVERGE  arcs;  and 
the  leftmost  arcs  are  ROTATION  arcs.  Goal  penalty  "costs"  must  be 
specified  for  all  the  ROTATION  arcs  and  GOAL  arcs.  We  have  chosen  to 
assign  a  penalty  of  2  units  for  rotation  at  the  right  time  to  the  wrong 
type  duty;  a  penalty  of  2  units  for  early  rotation  to  the  right  type 
duty;  and  a  penalty  of  4  units  for  late  rotation  or  early  rotation  to 
the  wrong  type  duty.  We  have  also  assigned  a  penalty  of  2  units  per 
person  for  under-  or  over-staffing  type  duty  1,  and  a  penalty  of  3  for 
under-  or  overstaffing  type  duty  2  or  3.  [These  penalties,  we  may  note, 
need  only  be  provisional.  They  can  be  used  to  obtain  a  trial  solution 
from  which  we  can  decide  whether  or  not  to  change  these  penalties  to 
better  reflect  priorities  or  policy  preferences,  as  we  shall  Illustrate 
with  an  additional  example.]  A  summary  of  the  actual  numbers  appears  in 


Type  duty 
From 


Time  on  Tour  1 
To  1  2  3 


Time  on  Tour  2 
1  2  3 


Time  on  Tour  3 
1  2  3 


1 

2 

3 


0  4  2 
2  0  4 
4  2  0 


4  2  0 
2  0  4 
4  2  0 


4  2  0 
0  4  2 
2  0  4 


Table  A-4  Goal  Costs  for  Rotation  Arcs 


The  above  Table  A-4  gives  the  simplified  data  for  this  hypothetical 
example.  The  penalty  for  "rotating"  from  Type  duty  1  to  Type  duty  1,  or 
Type  duty  2  to  Type  duty  2  etc.  (i.e.,  not  rotating  at  all)  during  Time 
on  Tour  1  is  thus  0,  since  no  one  should  be  rotated  with  Time  on  tour  1. 
The  penalty  for  rotating  from  Type  duty  1  to  type  duty  3  during  Time  on 
tour  1  is  2  units,  since  it  is  desired  that  personnel  on  Type  duty  1 
rotate  to  Type  duty  3,  although  not  before  completing  one  full  period  on 
tour.  The  penalty  for  rotating  from  type  duty  1  to  type  duty  2  before 
completing  one  full  period  is  4  units,  since  this  violates  two  aspects 
of  rotation  policy,  rotating  too  soon,  and  to  the  wrong  type  duty.  The 
other  penalties  are  similarly  explained,  and  all  the  penalties  for  this 
example  are  summarized  in  the  Table. 

Having  set  the  costs  on  all  the  arcs  i-.i  this  example,  we  need  to 
initialize  the  network  with  starting  inventories  of  personnel,  which  we 
assume  are  as  in  Tables  A-5.  Thus,  for  instance,  as  noted  in  this 
Table,  we  have  10  persons  with  Type  Duty  1  and  time  on  tour  1;  40 
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persons  with  Type  Duty  2  and  time  on  tour  1;  etc. 


Type  duty 

1 

2 

3 

Time  on  Tour 

1 

10 

40  . 

15 

2 

10 

40 

15 

3 

10 

40 

15 

Table  A-5  Starting  Inventories 

While  we  have  spread  the  personnel  out  evenly  over  the  various  tour 
lengths  to  simplify  computation,  in  the  real  GPSSR  the  distribution  of 
personnel  by  tour  length  will  be  determined  by  the  actual  data. 

Finally,  in  order  to  complete  the  network,  we  need  to  state  our 
desired  staffing  levels.  As  summarized  in  Table  A-6,  we  have  chosen  as 
goals  35  persons  in  Type  Duty  1,  125  persons  in  Type  Duty  2,  and  50 
persons  in  Type  Duty  3..  These  numbers  were  chosen  to  reflect  possible 
situations  where  goals  might  exceed  available  inventories  and  where  it 
is  not  possible  to  conform  to  stated  rotation  policies  in  all  detail. 
This  will  help  to  show  how  GPSSR  could  be  used  to  resolve  such 
conf 1 icts . 


Type  duty 
Staffing  Goal 


1  2  3 

35  125  50 


Table  A-6  Goals  for  number  of  personnel  in  each  type  duty 
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Since  our  hypothetical  inventory  is  less  than  our  desired  staffing 
levels,  it  will  not  be  possible  to  meet  all  goals.  The  model  solution 
obtained  from  GPSSR  is  "as  close  as  possible,"  however,  with  the 
resulting  divergence  shown  in  Table  A-7. 


Type  duty 

1 

2 

3 

Staffing  Goal 

35 

125 

50 

Best  Possible  Schedule 

35 

110 

50 

Divergence 

0 

15 

0 

Table  A-7  Goals  vs  Best  Possible 


Using  the  dollar  PCS  (personnel  transfer)  cost  from  Table  A-3  we 
find  that  the  solution  displayed  in  Table  A-7  will  involve  a  dollar  cost 
of  5570,000.  In  addition,  30  people  (about  15%)  had  to  be  rotated  to 
the  second  choice  in  the  preferred  rotation  sequence.  No  one  is  rotated 
either  early  or  late  as  a  result  of  this  optimization.  Some  of  the 
inventory  of  Type  1  duty  were  already  overdue  for  rotation,  but  all 
these  were  rotated. 

A  complete  list  of  all  rotations  and  costs  is  given  in  Table  A-9, 
but  before  considering  this  solution,  let  us  first  consider  whether  the 
"goal  penalty  costs"  that  yielded  this  solution  are  the  appropriate  ones 
to  meet  imperatives  in  the  rotation  pattern.  That  is,  we  want  to 
discover  whether  other  alternatives  might  be  preferable,  and  for  this 

purpose  we  want  to  bring  some  of  the  alternatives  into  view  in  an 
explicit  manner.  If,  for  example,  it  were  imperative  that  a_Q  the 
staffing  goals  for  Type  Duty  2  are  to  be  attained,  then  the  goal  cost 


for  Type  duty  2  should  be  set  to  some  larger  number,  like  10.  Such  a 
rearrangement  of  goal  deviation  penalties  results  in  a  new  rotation 
pattern  as  is  summarized  in  Table  A-8. 


Type  duty 

1 

2 

3 

Staffing  Goal 

35 

125 

50 

Best  Possible  Schedule 

35 

125 

35 

Di vergence 

0 

0 

15 

Table  A-8  Goals  vs  Best  Possible,  Type  Duty  2  Pre-emptive 


Comparing  the  rotations  summarized  in  Tables  A-7  and  A-8,  we  see 
that  25  persons  (about  12%)  had  to  be  rotated  to  the  second  choice  in 
the  rotation  sequence,  and  10  persons  (about  5%)  had  to  be  rotated  late. 
The  dollar  cost,  on  the  other  hand,  was  reduced  from  $570,000  to 
$405,000  in  going  from  Table  A-7  to  A-8,  since  the  staffing  goals  for 
Type  2  personnel  were  met  by  holding  back  personnel  due  to  be  rotated. 

To  conclude  this  example,  we  show  ,  in  Table  A-9,  the  complete 
solution  for  the  case  illustrated  in  Table  A-7,  with  all  the  rotations 
and  costs  for  each  type  of  duty  and  length  of  service.  Looking  at  the 
first  row  of  the  table,  the  FROM  NODE  is  characterized  by  Type  Duty  and 
time  on  tour,  as  is  the  TO  NODE.  In  other  words,  the  first  row  of  data 
is  from  Type  Outy  1 — tour  length  1,  to  Type  Duty  1 — tour  length  2. 

Since  this  is  the  prefered  transition,  the  GOAL  COST— i.e.,  the  penalty 
for  deviation  from  this  goal  —  is  0;  and,  since  no  PCS  move  is  involved 
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in  remaining  in  Type  Duty  1,  DOLLAR  COST  is  also  $0.  The  ARC  FLOW  is 
10,  indicating  that  all  10  persons  starting  in  Type  1  Duty  with  less 
than  1  full  period  service  (time  on  tour  1)  were  transfered.  ARC 
COST(G),  as  explained  earlier,  is  the  total  penalty,  i.e.  the  product  of 
GOAL  COST  and  ARC  FLOW.  In  this  case,  ARC  COST  is  0  since  no  penalty 
was  incurred.  ARC  C0ST($)  is  the  dollar  equivalent  of  ARC  COST(G),  and 
is  also  $0,  since  no  PCS  cost  was  incurred  on  this  arc.  Finally,  the 
column  labelled  MARG  COST(G)  indicates  the  rate  of  increase  of  goal  cost 
per  unit  increase  of  flow  (of  personnel)  on  an  arc  (for  small 
increases).  On  our  first  row,  the  absence  of  any  MARG  COST  indicates 
that  this  row  is  part  of  the  solution,  and  no  additional  cost  would  be 
incurred  by  adding  any  personnel  flow  on  this  arc,  i.e.,  bringing  this 
arc  into  the  solution,  since  it  already  is  in  the  solution.  The  values 
of  6  which  do  appear  in  later  rows  indicate  that  the  goal  cost  or 
penal ty  increase  per  additional  person  on  such  a  rotation,  should  such  a 
rotation  be  allowed  as  part  of  the  solution,  is  6.  This,  in  turn, 
implies  that  a  change  of  6  units  of  goal  costs  would  be  necessary  before 
this  arc  could  be  considered  for  inclusion  in  the  final  solution — i.e. 
before  any  personnel  would  be  considered  for  this  rotation. 
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FROM 

TO 

GOAL 

DOLLAR 

ARC 

ARC 

ARC 

MARG 

NODE 

NODE 

COST 

COST 

FLOW 

C0ST(G) 

C0ST($) 

C0ST(G) 

Type  Time  Type  Time 
Duty  Tour  Duty  Tour 

1 

1 

1 

2 

0 

0 

10 

0 

0 

1 

1 

2 

1 

4 

5000 

0 

0 

0 

1 

1 

3 

1 

2 

20000 

0 

0 

0 

1 

2 

1 

3 

4 

0 

0 

0 

0 

6 

1 

2 

2 

1 

2 

5000 

5 

10* 

25000 

1 

2 

3 

1 

0 

20000 

5 

0 

100000 

1 

3 

1 

3 

4 

0 

0 

0 

0 

6 

1 

3 

2 

1 

2 

5000 

10 

20* 

50000 

1 

3 

3 

1 

0 

20000 

0 

0 

0 

2 

1 

1 

1 

2 

5000 

0 

0 

0 

6 

2 

1 

2 

2 

0 

0 

40 

0 

0 

2 

1 

3 

1 

4 

9000 

0 

0 

0 

6 

2 

2 

1 

1 

2 

5000 

0 

0 

0 

6 

2 

2 

2 

3 

0 

0 

40 

0 

0 

2 

2 

3 

1 

4 

9000 

0 

0 

0 

6 

2 

3 

1 

1 

0 

5000 

25 

0 

125000 

2 

3 

2 

3 

4 

0 

0 

0 

0 

2 

3 

3 

1 

2 

9000 

15 

30* 

135000 

3 

1 

1 

1 

4 

20000 

0 

0 

0 

6 

3 

1 

2 

1 

O 

c. 

9000 

0 

0 

0 

3 

1 

3 

2 

0 

0 

15 

0 

0 

3 

2 

1 

2 

4 

20000 

0 

0 

0 

6 

3 

2 

2 

1 

2 

9000 

0 

0 

0 

3 

2 

3 

3 

0 

0 

15 

0 

0 

3 

3 

1 

1 

2 

20000 

0 

0 

0 

6 

3 

3 

2 

1 

0 

9000 

15 

0 

135000 

3 

3 

3 

3 

4 

0 

0 

0 

0 

6 

Total s 

195 

60 

570000 

*  indicates  a  penalized  rotation 

Table  A-9.  Actual  Flows  in  the  Optimized  Network 


A22 


OPTIMIZATION  AND  MONITORING  PROCEDURES 


The  operating  procedures  for  GPSSR  can  be  divided  into  two  categories, 
depending  upon  whether  an  optimization  function  (using  the  analytical 
goal -programming  network  model)  or  a  monitoring  function  is  desired.  The 
monitoring  function  can  also  be  considered  a  sub-function  of  the  optimiza¬ 
tion  function  as  well  as  an  (important)  function  in  its  own  right.  The 
procedures  presented  in  the  next  two  sections  are  those  necessary  for  exe¬ 
cuting  the  GPSSR  sequence  of  programs.  As  mentioned  earlier,  plans  are  in 
progress  for  more  "intelligent"  user-friendly  procedures  which  would 
require  minimal  knowledge  of  computer-related  concepts. 


OPTIMIZATION  PROCEDURES 

During  this  developmental  stage  of  GPSSR,  with  only  a  few  represen¬ 
tative  DCs  being  considered,  it  is  possible  to  maintain  these  DCs  data  on 
a  disk  file,  and  thus  "on-line."  The  production  version  of  GPSSR,  howev¬ 
er,  must  be  able  to  access  data  for  any  of  the  entire  set  of  DCs.  For  this 
reason,  the  data  for  all  DCs  will  reside  on  magnetic  tape.  Thus,  the  ini¬ 
tial  step  of  GPSSR' s  execution  must  be  one  of  reading  the  specified  DCs 
data  from  magnetic  tape  onto  disk,  creating  an  "on-line"  environment  for 
that  DC.  The  five  modules  of  GPSSR  must  be  executed  in  the  order  listed, 
thereby  restricting  the  user  to  issue  commands  in  a  specified  order,  as 
shown  below: 


MODULE  "Extraction  of  Data" 


EXTRACT  nn  :  Mounts  tape  of  DC  nn,  and  extracts  relevant  fields  on  to 
a  disk  file,  producing  an  "on-line"  environment  for  DC  nn. 


MODULE  "Calculation  of  Transition  Rates" 

TRANSIT  nn  :  Computes  transition  rates  for  DC  nn. 
GPSSMO  nn  :  Smooths  transition  rates  for  DC  nn. 


MODULE  "Network  Generator" 

GPSSRO  nn  :  Performs  the  advanced  start  for  the  network  opti¬ 
mizer  of  DC  nn.  Computes  staffing  levels  based  on  his¬ 
torical  and  user  defined  rates  (before  flexibility 
introduced) . 


MODULE  "Network  Optimizer" 

VICNET  nn  :  Computes  solution  to  sea/shore  rotation  problem 
for  DC  nn.  i.e.  minimizes  deviations  from  goals  while 
minimizing  actual  dollar  costs. 


MODULE  "Report  Generator" 

REP0RT1  nn  :  Produces  a  summary  report  for  DC  nn. 


Note  that  if  the  user  wishes  to  provide  input  in  the  form  of 
smoothing  or  flexibility  parameters  or  "overriding"  transition 
rates,  he  or  she  must  do  so  by  editing  (creating  or  modifying)  a 
NAMELIST  file  prior  to  the  execution  of  the  appropriate  task.  The 
format  of  the  NAMELIST  file  is  as  follows: 


SNAME  ALPHA=. 1 , 

DELTA=. 1 , 
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MONITORING  PROCEDURES 


A  critical  role  of  GPSSR  is  one  of  monitoring  past  and  present 
consequences  of  personnel  management.  Using  EMR  data,  the  system 
can  generate  a  variety  of  descriptive  statistics  and  display  them 
in  formats  which  are  meaningful  for  managers.  This  important  mon¬ 
itoring  function  can  be  achieved  by  means  of  the  three  "plotting" 
tasks,  directed  by  the  commands  PLOTSS,  PLOTTT,  and  PLOTSM,  as 
fol 1 ows : 

PLOTTSS  nn  :  Plots  personnel  at  sea  versus  personnel  on  shore  for 
DC  nn . 

PLOTTT  nn  :  Plots  time  on  tour  for  DC  nn. 

PLOTSM  nn  :  Provides  plots  of  smoothed  data  for  DC  nn. 

In  order  to  monitor  the  community,  the  optimizer  need  not  be  invoked,  nor 
must  the  monitoring  be  used  when  optimizing;  however,  the  extraction  mod¬ 
ules  discussed  in  the  previous  section  must  be  called  before  the  plotting 
routines.  The  next  section  shows  the  required  processing  order  of  the 
GPSSR  functions  and  modules. 


APPENDIX  B 


EXPONENTIAL  SMOOTHING 


Given  a  series  of  historical  rates,  which  may  be  trendy  and  noisy,  a 
common  method  for  estimating  a  "true"  current  rate  is  exponential  smooth¬ 
ing.  Proceeding  in  the  manner  of  an  exponential  function,  this  technique 
weights  current  data  more  heavily  than  the  earlier  data.  The  user  selects 
a  parameter,  a,  which  determines  how  much  additional  weight  should  be 
given  tc  the  current  year.  Choosing  a  =  0  gives  equal  weight  to  all  years. 
This  is  equivalent  to  taking  the  mean  of  the  time  series  as  the  estimate 
for  the  current  value.  Conversely,  choosing  a  =  1  uses  only  the  current 
year  as  the  estimate  of  the  "true"  value  for  the  series. 

In  our  case  the  time  series  consists  of  historical  rates  of 
promotion,  loss,  etc.,  for  each  of  the  past  4  years.  We  need  an  estimate 
of  the  rates  for  the  next  4  or  5  years.  The  rates  for  even  a  stable  commu¬ 
nity  tend  to  oscillate  somewhat,  and,  when  the  oscillation  is  not  too 
great,  as  is  true  for  most  stable  communities,  the  mean  rates  would  be 
most  appropriate.  For  an  expanding  community,  however,  especially  a 
community  which  has  started  expanding  less  than  4  years  ago,  the  mean 
rates  are  less  appropriate  than  a  weighted  series  in  which  the  most  recent 
year  is  given  greater  weight. 

The  exponential  smoothing  formula  is  not  generally  given  in  closed 
form,  but  is  usually  given  recursively.  If  R(n)  is  the  rate  for  year  n, 
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and  S(n-l)  is  the  smoothed  rate  for  year  n-1,  then  the  formula  for  S(n)  is 
given  by 

S(n)  =  a  *  R(n)  +  (  1  -  a  )  *  S(n-l) 

For  the  first  year 


S(l)  =  R(  1) 

In  our  program,  we  actually  take  a  modified  S(n),  S'(n),  where 

S'(n)  =  a  *  S(n)  +  (1  -  a  )  *  M 
where  M  is  the  mean  rate  for  the  entire  series.  1 


As  an  example,  suppose  the  data  for  5  years  are  10,8,11,9,12.  This 
represents  a  series  8.5,9,9.5,10,10.5  with  "noise"  of  1.5,— 1,1.5,— 1,1.5 
"added".  Choosing  a  =  0  gives  an  estimate  of  S  =  10,  which  is  too  low  an 
estimate  for  the  current  average  value  of  the  series.  Choosing  a  =  1 
gives  an  estimate  of  S  =  12,  which  is  too  high.  The  next  value  of  the 
series  will  actually  be  10,  but  with  a  "true"  value  of  11.  An  a  of  .3  gives 

the  estimate  S  =  10.10,  while  an  a  of  .1  gives  an  estimate  of  10.005.  ,  Both 
of  these  are  a  better  estimate  than  the  mean  for  the  "true"  value  of  the 

series . 

1  The  user  is  allowed  to  vary  a  for  each  rate.  There  is  no  compelling 
unique  choice  of  a  ,  beyond  the  obvious  observations  that,  if  the  communi¬ 
ty  is  stable,  a  =  0  is  the  correct  choice,  while  if  data  earlier  than  the 
current  year  is  irrelevant,  then  a  =  1  should  be  used.  Many  textbooks 
recommend  choosing  a  between  .01  and  .3.  However,  we  have  run  a  number  of 
tests  with  different  choices  of  a  ,  but  without  any  decisive  results.  The 
user  not  skilled  in  time  series  analysis  is  advised  to  use  the  default 
values  for  a . 
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